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Abstract 

This is the first of a series of papers extending a 1+3 covariant and and gauge invariant 
treatment of kinetic theory in curved space-times to a treatment of Cosmic Background Radi- 
ation (CBR) temperature anisotropics arising from inhomogencities in the early universe. This 
paper deals with algebraic issues, both generically and in the context of models linearised about 
Robertson- Walker geometries. 

The approach represents radiation anisotropics by Projected Symmetric and Trace-Free ten- 
sors. The Angular correlation functions for the mode coefficients are found in terms of these 
quantities, following the Wilson-Silk approach, but derived and dealt with in 1+3 covariant and 
gauge invariant (CGI) form. The covariant multipole and mode-expanded angular correlation 
functions are related to the usual treatments in the literature. The CGI mode expansion is re- 
lated to the coordinate approach by linking the Legendre functions to the Projected Symmetric 
Trace-free representation, using a covariant addition theorem for the tensors to generate the 
Legendre Polynomial recursion relation. 

This paper lays the foundation for further papers in the series, which use this formalism in 
a CGI approach to developing solutions of the Boltzmann and Liouville equations for the CBR 
before and after decoupling, thus providing a unified CGI derivation of the variety of approaches 
to CBR anisotropies in the current literature. 
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1 Introduction 

Ellis, Treciokas and Matravers (ETM) ||, |j introduced a 1 + 3 covariant kinetic theory formalism 
in which an irreducible representation of the rotation group based on Projected Symmetric and 
Trace-Free (PSTF) tensors orthogonal to a physically definved 4-vrelocity u a gives a covariant 
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representation of the Cosmic Background Radiation (CBR) anisotropies, which is gauge-invariant 
when the geometry is an almost-Robertson Walker (RW) geometry. This 1 + 3 Covariant and 



Gauge-Invariant (CGI) formalism has been used in a previous series of papers [||, [jl], 11, [0| to 
look at the local generation of CBR anisotropies by matter and spacetime inhomogeneities and 
anisotropies in an almost-Friedmann Lemaitre (FL) universe model []. By contrast, the present 
series of papers^ uses this formalism to investigate CBR anisotropies in the non-local context of 
emission of radiation near the surface of last scattering in the early universe and its reception here 
and now (the Sachs- Wolfe (SW) effect and its further developments). 

There is of course a vast literature investigating these anisotropies both from a photon view- 
point, developing further the methods of the original Sachs- Wolfe paper 2£], and from a kinetic 
theory viewpoint, so it is useful to comment on why the CGI philosophy and programme make 
the present series of papers worthwhile. Rather than beginning with a background described in 
particular coordinates and perturbing away from this background, this approach centres on 1 + 3 
covariantly defined geometric quantities, and develops exact nonlinear equations for their evolution. 
These equations are then systematically linearised about a Friedmann-Lemaitre (FL) background 
universe with a Robertson- Walker (RW) geometry resulting in description by gauge-invariant vari- 
ables and equations [16]. Because the definitions and equations used are coordinate-independent, 
one can adopt any suitable coordinate or tetrad system to specialise the tensor equations to specific 
circumstances when carrying out detailed calculations; a harmonic or mode analysis can be carried 
out at that stage, if desired. 

This approach is geometrically transparent (see |30|, j?l]]) because of the CGI variable definitions 
used. In contrast to the various gauge-dependent approaches to perturbations in cosmology, the 
differential equations used are of just the order that is needed to describe the true physical degrees of 
freedom, so no non-physical gauge modes occur. When a harmonic decomposition is introduced in 
the case of linear perturbations, the CGI variables used here provide a description that is equivalent 
to that obtained by approaches based Bardeen's GI variables [33], see [34], but they do not imply 
linearisation of the equations from the outset, as occurs in that formalism. 

Thus the benefit of the present formalism is precisely its 1 + 3 covariant and gauge invariant 
nature, together with the fact that we are able to write down the exact non-linear equations 
governing the growth of structure and the propagation of the radiation, and then linearise them in a 
transparent way in an almost-RW situation. This means it can be extended to non-linear analyses in 
a straightforward way [14], which will be essential in developing the theory of finer CBR anisotropy 
structure as reliable small-angle observations become available. Achievements of the CGI approach 
with respect to the CMB are the almost-EGS Theorem |J, related model-independent limits on 
inhomogeneity and isotropy |ll|, [l^, |13| , and derivation of exact anisotropic solutions of the 
Liouville equation in a RW geometry ( , see also ) . 

This paper, Part I, deals with algebraic issues, developing further the formalism of ETM: namely 
an irreducible representation of radiation anisotropies based on PSTF tensors ||, pj. The paper 
considers this irreducible representation and its relation to observable quantities, both generically 
and in the context of models linearised about RW geometries JO]]. In section 2 and 3, the underlying 
1 + 3 decomposition is outlined and the basic CGI harmonic formalism for anisotropies developed. 
In section 4, the angular correlation functions are constructed from CGI variables, assuming that 
the multipole coefficients are generated by superpositions of homogeneous and isotropic Gaussian 
random fields. The multipole expansion is discussed in detail, extending the results of ETM, 

1 Here 'Robertson- Walker' refers to the geometry, whatever the field equations; 'Friedman-Lemaitre' assumes that 
the Einstein gravitational field equations with a perfect fluid matter source are imposed on such a geometry. 

2 And a series by Challinor and Lasenby po[ |, |3l| which are similar in method and intent but different in detail 
and focus. 
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giving the construction of the multipole coefficient mean-square and developing its link to the 
angular correlation function. In section 5, the mode coefficients are found following the Wilson- 
Silk approach, but derived and dealt with in the CGI form; the covariant and gauge invariant 
multipole and mode expanded angular correlation functions are related to the usual treatments 
used in the literature [17, [20, E^. 



In this discussion, the CGI mode expansion is related 
to the coordinate approach by linking the Legendre Tensors to the PSTF representation, using a 
covariant addition theorem to generate the Legendre Polynomial recursion relation. The key result 
is the construction of the angular correlation functions in the CGI variables, and their link to the 
(non-local) GI Mode functions [20]. 

The following papers in the series look at the Boltzmann equation and multipole divergence 
relations, solution of the resulting mode equations, and relation of the kinetic theory approach to the 
photon based formalism of the original Sachs- Wolfe paper. Exact non-linear equations are obtained 
and then linearised, allowing a transparent linearisation process from the non-linear equations that 
is free from ambiguities and gauge modes. 



2 Temperature anisotropies 

A radiation temperature measurement is associated with an antenna temperature, T(x*,e a ), mea- 
sured by an observer moving with 4-velocity u a at position x % in a direction e a on the unit sphere 
(e a e a = 1, e a u a = 0). We assume u a can be uniquely define in the cosmological situation, corre- 
sponding to the motion of 'fundamental observers' in cosmology |37j The direction e a can be 
given in terms of an orthonormal tetrad framed], for example by : 

e a (9,(j)) = (0, sin 6 sin </>, sin 6 cos 4>, cosO) . (1) 

The temperature T(x t ,e a ) can be unambiguously decomposed into the all-sky average bolometric 
temperature T(x l ) at position x l , given by 

T(x*) = ±- / T(x\e a )dQ. (2) 

where f2 is the solid angle on the sky, and the anisotropic temperature perturbation 5T(x l ,e a ) (the 
difference from the average over the unit sphere surrounding x % ]l0|])> can be defined : 

T{x\e a ) =T{x l ) +5T{x\e a ). (3) 

From the Stefan-Boltzmann law it follows - if the radiation is almost black-body, which we assume 
- that the radiation energy density is given in terms of the average bolometric temperature by 
Pr{ x1 ) = rT 4 (x*). (r is Stefan-Boltzmann constant). Both the quantities T(x l ) and 5T(x l ,e a ) are 
CGI, for T(x l ) is defined in a physically unique frame in the real universe (because u a is assumed to 
be uniquely defined), and 5T(x l , e a ) vanishes in any background without temperature anisotropies. 
We can define the fractional temperature variation t(x 1 , e a ) by [[H]] 

3 In the early universe, when matter and radiation average velocities differ, there may be several competing 
possibilities for covariant definition of u a ; however once a choice has been made between these possibilities, this 
vector field is uniquely defined. 

4 The description of such a tetrad frame is briefly discussed in last appendix 

5 Note that this is not the same as the background temperature, for that quantity varies only with cosmic time t, 
whereas the true isotropic component of the temperature varies with spatial position as well as time. 
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and take a covariant (angular) harmonic expansion of this, 

oo 

r(x\e a ) =Y,Taiaza 3 ...a t (x i )e^e a *e a K..e a ^e a i e^t^. (5) 

1=1 £>1 

We introduce the shorthand notation using the compound index An = a\ai...ai. Here r aia2a3 ... a; (x l ) 
are trace-free symmetric tensors orthogonal to u a : 

T A e = T(A e ) i TA e abh ab = , TA ia u a (6) 

Round bracketts "(••)" denote the symmetric part of a set of indices, angle bracketts "(..)" the 
(orthogonally-) Projected Symmetric Trace-Free (PSTF) part of the indices : ta z = T (A t )- 

Because of (||), this expansion is entirely equivalent to a more usual expansion in terms of 
spherical harmonics: 

oo 

r(x i ,e a ) = Y, A T(^)yr(0,4>) (7) 
1=1 

(see |U for details), but is more closely related to a tensor description, and so results in more 
transparent relations to physical quantities. 

We wish to measure the temperature in two different directions to find the temperature differ- 
ence associated with the directions e a and e' a such that (using ||) : 

AT(x i ; e a , e' a ) = T(x\ e a ) - T(x\ e' a ) => AT(x'; e a , e' a ) = 5T(x\ e a ) - 5T{x\ e' a ). (8) 

It follows from (§), (|) and (|) that AT(x^; e a , e' a ) = T(x i )Y, e .r Al (e A * - e' A A , where AT/T 
represents the real fractional temperature difference on the current sky. Due to the CGI nature of 
T(x l ) we may relate this directly to the real temperature perturbations (no background model is 
involved in these definitions). 

The relation between the two directions e a and e' a at x % is characterised by 

e a e' a = cos(/3) =: X. (9) 



l a 



i.e. they are an angular distance (5 apart [J If analogous to (|l|) we write e a (9',(J)') = e 
(C^sin^'sin^sin^'cos^cos^'), then it can be shown from (H|) that 

cos (3 = sin 6 sin 9' cos{4> — <p') + cos cos 9 ' . (10) 

In later applications, it is important to relate the different terms of the harmonic expansion to 
angular scales in ths sky. A useful approximation is I ~ a, where 6 is in radians. 



2.1 Covariant and gauge invariant angular correlation function 

The two-point correlations are an indication of the fraction of temperature measurements, T(x l , e a ), 
that are the same for a given angular separation. This corresponds to the correlation between 
5T(x l ,e a ) and 8T(x l ,e' a ) or equivalently between r(x l ,e a ) and r(x l ,e /a ), given by the angular 
position correlation function 

C(e a ,e' a ) = (r(x\e a ),r(x\e ,a )), (11) 

6 It should be pointed out that e' a is distinct from e a , the first denotes a direction vector different from e a in a 
given tetrad frame, while the second means the same direction vector in a different tetrad frame. 
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where the angular brackets representing an angular average over the complete sky. Note this is a 
function in the sky. If we write t(x 1 , e a ) and t(x j , e' a ) in terms of the angular harmonic expansion 
(H), we can also define correlation functions C\ for the anisotropy coefficients TA l (e a ),TA l (e' a ) by 

C e (e a , e' a ) = (2£ + l)' 1 A, (r Al (x\ e a )r A \x\ e a )) . (12) 

Here the right-hand side term in brackets is the all-sky mean-square value of the l—ih temperature 
coefficient TA e (x l , e a ), and the coefficient A^ is defined in ( |119| ). The numerical factor (21 + 1) A^ 
is included in order to agree with definitions normally used in the literature (see later). This can 
be thought of as the momentum space version of ( pd| ) , as we have taken an angular fourier series of 
the quantities in that equation; it says, for each choice of e a ,e' a , how much power there is in that 
expression for that angular separation as contributed by a particular £-ih valued multipole moment 
on average. 

2.2 The Central-Limit Theorem 

We consider an ensemble of temperature anisotropies, where a sequence of repeated trials is replaced 
by a complete ensemble of outcomes. The temperature anisotropy Tj± t found in a given member of 
the ensemble is a realization of the statistical process represented by the ensemble. The physically 
measured anisotropy is taken to be one such realization. The variance of the ensemble, for example 
(rA e T Ae ), is in principle found by averaging over a sufficiently large number of experiments, where 
we assume the results will approach the true ensemble variance - this is the assumption of ergodicity. 

On Fourier transforming, we make the assumption that to a good approximation the phases 
of the various multipole moments are uncorrelated and random. This corresponds to treating the 
anisotropies as a form of random noise. The random phase assumption has a useful consequence: 
that the sum of a large number of independent random variables will tend to be normally dis- 
tributed. By the central-limit theorem this is true for all quantities that are derived from linear 
sums over waves. The end result is that one ends up with a Gaussian Random Field (GRF) which 
is fully characterized by a power spectrum. The central limit theorem holds as long as there exists 
a finite second moment, i.e. a finite variance. 

We will assume that the angular variance (t^t^) is independent of position; this is the assump- 
tion of statistical homogeneity. Its plausibility lies in the underlying use of the weak Copernican 
assumption. Additionally it is convenient to assume that the power spectrum will have no direc- 
tional dependence, thus it will be isotropic: P(k a ) = P(\k a \). Together these imply the statistical 
distribution respects the symmetries of the RW background geometry. 

One needs to be careful in using the central limit theory to motivate Gaussian random field, 
particularly in the presence of nonlinearity, which could result in the elements of the ensemble no 
longer being independent. While the assumptions of primordial homogeneous and isotropic GRF's 
is plausible, because the perturbations are made up of a sufficiently large number of independent 
random variables, the key point is to realize that these are assumptions that should be tested if 
possible. The simplest test of weak non-Gaussianity is looking for a three-point angular or spatial 
correlation, for the Gaussian assumption ensures that all the odd higher moments are zero and that 
the even ones can be expressed in terms of the variance alone. If the primordial perturbations are 
made up of GRF's, then non-Gaussianity of the CBR anisotropy spectrum should arise primarily 
from foreground contamination due to local physical processes. If the non-Gaussian effects due to 
these later physical processes or evolutionary effects are small enough, one can attempt to determine 
a cosmological primordial signature. 
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2.3 Gaussian perturbations 

A general Gaussian perturbation [27], r(x*,e a ), will be a superposition of functions, T Ae , i.e. (||) is 
satisfied, where the probability, P, of finding a particular valued temperature coefficient is given by 
(a* = (t Ai t a *)) : 

F (^) = T= ex p{ Z ^fi^l- ( 13 ) 

Note that r Ae is both the amplitude of the £-th component, and determines the probability of 
that amplitude. The probability of a temperature perturbation, r, is given by the sum of the 
Gaussian probability distributions (|l^) weighting the various angular scales, given by I, of the 
general perturbation (||). 

Considering isotropic and homogeneous gaussian random fields, the angular position correla- 
tion function C(e a ,e >a ) is a function only of the angular separation (3 of the two temperature 



measurements. We then write (11) as 

(r • t') = C((3) = W(X) , (14) 

where the expression on the left is shorthand for (r(x l , e a ), t(x 1 , e fa ))p , the 2-point angular cor- 
relation function for a given angular separation (5 between the on-sky temperature measurements, 
and X = e a e' a = cos (3. This expression is now independent of position in the sky. Gaussian Fields 
are completely specified by the angular power spectrum coefficients Ce (|i~2]), which are now just 
constants, because I is uniquely related to /?, so the power spectrum is a function of the modulus 
of the wavenumber only. One thus expects the temperature perturbations in this case to be fully 
specified by the mean squares, (rA e T Ali ^j, when ^ is substituted in (|il|). Equivalently they are 
uniquely determined by the angular Fourier transform of the 2-point angular correlation function . 



3 Multipole expansions 

In this section we examine the anisotropy properties of radiation described in terms of the covariant 
multipole formalism (||), which is equivalent to the usual angular harmonic formalism but much 
more directly related to space-time tensors. Note that the relations in this section hold at any point 
in the space-time, and in particular at the event R ('here and now') where observations take place. 
Here we consider the PSTF part of, e Al ; some useful properties of e Al are listed in appendix B. 

3.1 The PSTF part of e Ae 

Because the coefficients in (||) are symmetric and trace-free, the important directional quantities 
defined by directions e a at a position x % are the PSTF quantities 

A e = e (A e ) = e (a le a 2e a 3 ^ e a e _ le a e )^ ^ 

for clearly 

oo oo 

T(x\e a ) = Y j r Al ^ i )e Al {0^) = £ r Ae {x 1 )O a ^6, 0). (16) 

£=l 1=1 

Indeed the standard spherical harmonic properties are contained in these quantities. 

Now the symmetric trace-free (STF) part of a 3-tensor is given in general by Pirani ||], 

[^/ 2 ] (—l) n £\(2£ — In — 1)!! 

\. F A e ] =J2 Blnh (^2-- h a2n-ia2n F a2n+i:.a e ) with B tn = (£_ 2n )!(2^-l)!!(2n)!!' ^ 
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Here [£/2] means the largest integer part less than or equal to £/2. The following definitions have 
also been used : l\ = £{£ - 1)(£ - 2)(£ - 3)...(1), and £\\ = £(£ - 2)(£ -4){£- 6). ..(2 or 1). 
The PSTF part of a tensor, 



T {ab) = [T ab ] PSTF =[[T ab r"\' , (18) 

can be constructed recursively from a vector basis following EMT and a little algebra. 
We take the PSTF part of |, |, |, § to find 

[1/2] 

Q A t = e {A e ) = B ek h (A2k e A '- 2k \ (19) 

where h^ A2k e Al ' 2k) = /i (ai " 2 ...h a2k ~^ a2k e a2k + 1 e CL2k + 2 ...e ai) and B ik are given by (|l7|). 

From || we can now construct recursion relations that play a key role later on. First, 

Q(A t a t+1 ) = q(Ai at+i) — e a h ad O d ^- 1 /» a <°'+ l) (20) 

2^ + 1 v ; 



iSTF 



P 



From (p0|), (|l5|), and using 



e « a4 ' = (27rT) ""' < 21 » 



it can then be shown that 



p2 

U ~ {2£ + \){2£-l) n U [ZZ) 

relates the {£ + !) — th term to the £ — th term and the (£ — 1) — th term. 

The orthogonality, addition theorem and double integral relations of Al are listed in appendix 
B. Using the orthogonality relations we obtain the inversion of the harmonic expansion : 



^VJ r. 

T(x\e a ) = Y / TA i (x i )0 Ae O r Ae (x i ) = Aj 1 / dttO At T{x i ,e a ). (23) 
1=0 MlT 

The polynomial L t = Ai O' Ae = ES^" 2 " 1 - ( see S)) is the natural polynomial that 
arises in the PSTF tensor approach (equivalent to the Legendre polynomials, see below), where the 
coefficients B\ m are defined by (|l7|). It follows from this that 

Li(l) = A *0 Al = J2 B lm =: e with p t = ( = j^-^- (24) 



m=0 

The fits satisfy the recursive relations 



(^ + l) Wi ' (2£-l)^ _i ' (2£ + l)(2£-l) 

Any function W(X) can be expanded in terms of the polynomials L(_(X) [] and then upon 
combining ( p.24j ) and expansion in terms of L#(X) to find the expansion in terms : 

oo oo 

W(X) = £C^pQ and W(X) = Y / CnO A "0' An . (26) 

£=1 n=l 

When W(X) is the angular correlation function, the Ct are the corresponding angular power 
spectrum coefficients (see below). 

7 We include the dipole, I = 1, however the monopole is dropped as we will only be considering the expansion of 
the CGI perturbations ^ where the isotropic part has been factored out according to (^). We must beware that 
this does not cause problems later by omitting the spatial gradients of the isotropic term. 
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3.1.1 Relationship to Legendre polynomials 



A Legendre polynomial Pe(X) is given by renormalising the polynomials L^{X) denned in (|124| ) so 
that P e (l) = 1. By (M), this implies 



P l (X) = ((3 l y l L l {X) =>P t (l) = l; 



consequently from ([[2 



A '0' Ae = faP e (X), 
where (5i are given by (||). It follows from (||) that 



[t/2] 

Pt{X) = ]T A im X 
m=0 



with A tk - (-^W-W 



2 e k\(£ - k)\(£ - 2k)\ 



are related to the Ba~ in (123) by 



B 



Ik 



(2i)\ 



A 



(27) 
(28) 

(29) 
(30) 



Any function W(X) can be expanded in terms of both sets of polynomials - see ( |26|) and the 
corresponding expression 



OO cyn . -1 



(31) 



These two expansions can then be related as follows: equating ( p6| ) and (31), and using (pi) and 
(IH), gives 



oo [//2] 

E E CiBi m X 
l=\ m=0 



^ [ ^ ] + 1; 



EE 



4tt 



CiAi m X 



' CgBim = - — — ^-CiAim, (32) 



£=1 m=0 

from ( p0[ ) this gives the relation between the expansion coefficients : 

"(2tf + l)(2£)!" 



4tt 



4^2^!) 5 



(33) 



3.2 The mean square of PSTF coefficients : {F Al F At ) 

It is known as before, from evaluating J dilf 2 and constructing the orthogonality conditions on 
A( , that inversion 



„ oo 

F Ae = Af <MO A J(x\e a ) o f(x\e a ) = Y J FA l O Al , 



can be constructed. From this we can build 



F Ai F B ™ = AJ 1 A~ 1 [ cMO A J(x\e a ) f dtl'O lB ™f(x\ e' a " 

Jin J4n 



to find 



F Ap F m = Aj A m 



dndtt'0 Ae O' Bm f{x\ e a )f(x\ e' a ) . 



(34) 



(35) 



(36) 
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(41) 



Taking the ensemble average [[??], 18] then gives 

(F Ae F B ™) = A7 1 A m 1 | Jdndn'o Ae o ,B -(f-f) (37) 

where (...) indicates an ensemble average over sufficiently many realizations of the angular correla- 
tion function. 

In order to evaluate this further we assume that the correlations between the function f(x l , e a ), 
i.e. , {f(x l ,e a )f(x l ,e' a )) are a function of the angular separation between the two directions only, 

(f.f')=W(e a e' a ) = W(X). (38) 

This is a consequence of the Gaussian assumption (|i~4|), which allows one to evaluate the angular 
correlation functions ( |IT|) in a straightforward way 0. With this assumption (37) becomes 

(F Ae F B ™) = AfA- 1 J J dn<m'0 Ai O' Bm W(X) . (39) 

Substituting ( p6| ) into (|39l), we find 

(F Al F Bm ) = A^ 1 A m 1 / / dndn'0 A ,0' B ™ £ C n O An O' An ■ (40) 

n=l 

Rearranging terms, 

(F At F B ™) = AfA' 1 £ C n | dfiO^O A »| | jf dn'o'^o'^ 
where the integrals can be evaluated using the orthogonality conditions on the Oa/s, ( |1T 

oo 

(F Al F Bm ) = AfA m 1 Y J C n {5" n A,h {Al) ^}{5 mn At 1 h^ {An ^ 

n=l 

oo 

= J2 d n5 £n 5 mn h {Ae) ^h^\ An) . (42) 

n=l 

Thus 

(i%F B ») = C e S em h {A / B ™\ to find (F Ae F Ae ) = Cth^ A ^ At \ (43) 
so using fo/^)^' (121) the mean-square is found to be 

(F Ai F A *) = C e (2e + l) (44) 

giving the angular power spectrum coefficients in terms of the ensemble-averages of the harmonic 
coefficients 0. If we use the Legendre expansion (|3l] ) instead of the covariant expansion coefficients 
fl2~6|), then from ( p4| ) and (33) the relation is 

(F Ae F Ae ) = (21 + 1)A^ 1 Q, (45) 

where are the usual Legendre angular power spectrum coefficients. 

8 Ideally one would prefer to evaluate the angular correlation function without using the Gaussian assumption, as 
this is one of the features one should test rather than assume. 

9 This corrects an error in removing a spurious factor of 3 £ = h^ e which follows from the orthogonality 

conditions in la] which are corrected here. 
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3.3 The CGI angular correlation function 



We can now gather the results above in terms of the application we have in mind, namely anisotropy 
of the CBR. Consider on-sky perturbations made of Gaussian Random Fields : The angular cor- 
relation function C{fi) = W(X) is given by (14); the angular power spectrum coefficients Ci are 
given by the mean-square of the I— th temperature coefficient through (f45|): 

(r Al T Ae ) = (2^ + l)A7 1 Q (46) 



where the constants are given by ( 119 ), These quantities are related by (pl|): 

W(X) 



Cm- ~Pm{X) 



m=l 



47T 



(47) 



where the P m {X) are given by (|27|) from (31) and rt2 



3.3.1 Cosmic variance 

The observations are in fact of af (which is 2~2m=-e l^ml 2 /^ 71 " in the usual notation). This is 
what is effectively found from experiments, such as the COBE-DMR experiment. This is a single 
realization of the angular power spectrum Cg. The finite sampling of events generated by random 
processes (in this case Gaussian random fields) leads to an intrinsic uncertainty in the variance even 
in perfect experiments - this is sample variance, or in the cosmological setting, cosmic variance. 
We are measuring a single realization of a process that is assumed to be random; there is an error 
associated with how we fit the single realization to the averaged angular power spectrum. 

The quantities (TA e T Ae ) represent the averaged (over the entire ensemble of possible C/s) angu- 
lar power spectrum, this is what one is in fact dealing with in the theory, as the reductions are done 
in terms of Gaussian Random Fields where the entire ensemble is considered rather than a single 
experimental realization. The af are a sum of the U + 1 Gaussian Random Variables a^ m , this is 
taken to be x 2 distributed with 2£ + 1 degrees of freedom. Each multipole has 2£ + 1 samples^. 

The key point here is that cosmic variance is proportional to l~ l l 2 and so is less significant 
for smaller angular scales than larger scales (as is popular wisdom), i.e. cosmic variance is not an 
issue on small scales. However on small (and perhaps intermediate scales) systematic errors could 
be underestimated. 

Physical process deviations and instrument noise are expected to dominate the small scales 
rather than non-Gaussian effects in the primordial perturbations, but on large scales the uncertainty 
due to cosmic variance would swamp out a non-Gaussian signature. It then seems plausible that 
on both large and small scales the assumption of Gaussian Perturbations is acceptable; however on 
intermediate scales this is not the case, on these scales the effects of comic variance would be small 
enough to allow a non-Gaussian signature to be apparent. 



4 Mode expansions 

We now consider spatial harmonic analysis of the angular coefficients discussed in the previous 
section. Note that the relations in this section hold in space-like surfaces, namely the background 
space-like surfaces in an almost-FL model. The application in the following sections will be to the 
projection into these spacelike surfaces of null cone coordinates associated with the propagation of 
the CBR down the null cone. 

10 The uncertainty in Ce as — ./ , 2 i +1 \ 
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Following the Wilson-Silk approach [jTj, 20, |2^| we consider the following CGI expansions. 
Eigenfunctions Q{x v ) are chosen to satisfy the Helmholtz equation 

D«D a Q = -k 2 phys Q (48) 

in the (background) space sections of the given space-time of interest, where the Q's are time- 
independent scalar functions with the physical wavenumber k p h ys {t) = k/a(t), the wave number 
k being independent of time Fl These define tensors Q Ae (k u ,x' 1 ) that are Projected, Symmetric, 
and Trace-Free, and in the case of scalar perturbations are chosen to be given by PSTF covariant 
derivatives of the eigenfunctions Q: 

Q Al = {-k phys )- l V {Al) Q (49) 

Using these we define functions of direction and position: 

G e [Q}(x»,e a ) = A 'Q Ae (50) 

with the At defined by (|l5|) . Here the Gg are call mode operators and the objects Gg[Q] are called 
mode functions 0. It follows that 

G t [Q] = (-k phys )- e O A 'D {Ae) Q (51) 

and we can expand a given function f(x l , e a ) in terms of these functions. In our case this serves as 
a way of harmonically analysing the coefficients t Ai (x 1 ) in (||) and ([u]): expanding the temperature 
anisotropy in terms of the mode functions, 

r(x\e a ) = Y / T, T ^k)G e [Q] (52) 

i k 

where the ^-summation is the angular harmonic expansion and the /j-summation the spatial har- 
monic expansion (in fact k will be a 3- vector because space is 3-dimensional, see below). Using the 
expansion (||) on the left and (^) on the right, 

Y,TA t {x l )0 Al = EE r K*^)(-Wr^D ( ^)Q (53) 

i Ik 

and so 

t a ,(x<) = J2 T dt,k)(-k phys )- e D {Ae) Q(xn (54) 

k 

which is the spatial harmonic expansion of the radiation anisotropy coefficients in terms of the 
symmetric, trace-free spatial derivatives of the harmonic function Q. The quantities Te(t,k) are 
the corresponding mode coefficients. Note that we have not as yet restricted the geometry of the 
Q's: they could be either spherical or plane-wave harmonics, for example. 
By successively applying the background 3-space Ricci identity, 



K 

V abAt Q - D baAi Q = + E - { Sb b n ha n a - firh^i) D Ai Q, (55) 



n=l 



n The function Q will be associated with a direction vector e" fe ) and wave vector k a = ke^ normal to the surfaces 
Q — const, see the following subsection. 

12 Note these are functions in phase space, not on M. 
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where At = a\...a n -\b n a n j r \...ai, i.e. , the sequence of £ indices with the n-th one replaced with a 
contraction. First, the curvature-modified Helmholtz equation is found [^] (using 133[ , 134 ,132) : 



D a D a Q 



-k^Q 



with 



\ (K£{£ 



+ 2)-t 



Second, we are able to construct the mode recursion relation (using |130| , |l3l| in |2 

Gi+i[Q] 



e a D a [G,[Q]] = +k phys 



(2l + l)(2£-l) 



K 



l))Gi-i[Q] 



(56) 



(57) 



The latter is the basis of the standard derivation of the linear-FRW mode heirarchy for scalar 
modes. The derivation of these are given in appendix C (we consider only scalar eigenfunctions). It 
will be seen, in paper II, that this relation can be used in place of the general divergence relations 
which allow the construction of generic multipole divergence equations 0j if one restricts oneself to 
constant curvature space-times. 

Given the recursion relation one can immediatly make the connection with the usual Legendre 
tensor treatment (139) this is shown in appendix C (139-141). 



4.1 | Ti | 2 in almost-FLRW universes 

We now relate the multipole mean-squares (ta^ 1 ^ of the ensemble average over the multipole 

moments with that of the mode coefficient mean-squares |t£(/c)| 2 . 

In order to carry this out we relate two separate spatial harmonic expansions (|54]) for the same 
function: the first is one associated with plane wave harmonics (Q k ), naturally used in describing 
structure existing at any time t, and the second, one associated with radial and multipole harmonics 
(O^H Rai), i-e. , a spherical expansion based at the point of observation, naturally arises when 
we project the null cone angular harmonics into a surface of constant time. These are both related 
to the Mode Function formulation which becomes useful in the non-flat constant curvature cases. 



4.1.1 Plane- waves and Mode functions 

Considering flat FRW universes, each set of eigenfunctions satisfy (fief). The temperature anisotropy 
(Q) can be expressed in terms of its plane-wave spatial Fourier transform : 

r(x\e a )=Y / r(k,t,e^,e a )Qfiaf (58) 

(For a more detailed treatment see appendix D). It can be shown that for the flat case, K = 0, 
that is (Hp hold in jj^) and (|l43;) to find from (|l45|) : 

V( Ae )Qflat = {-ikphysYO^Qfiat (59) 

where are the PSTF tensors associated with the direction e^ k y Then from ( f49| ) we find that : 

QflatAt = (-l) E of f Q flat , (60) 



along with (|54|) to get the temperature multipole : 



T Ae = (-lYJ2Mt,k)Q flat O^. (61) 

k 



13 This will be necessary in order to switch from the almost-FLRW multipole divergence equations in part II to the 
mode representation where we use ( ^l|) to construct D a D ( a QA,)- 
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4.1.2 Radial expansions and Mode functions 

Using the flat, K = 0, spherical eigenfunctions centred on a point x l , and with associated radial 
direction vector ei x \ The latter is the same as the (spherically symmetric) projection into the 
constant time surfaces of the tangent vector e a of the radial null geodesies, so we need not distinguish 
it from that vector. In this case the £-th harmonic is 



Ql(x a ,e a {x) )\ x = R Ae (r)0^, D a r = e a , D a e a 



(62) 



where e a = dx a /dr is the unit radial vector. Cartesian coordinates in space are given by r and e a 
through x l = re 1 . Defining the projection tensor 



Pab = Kb ~ e a e b pabu a = = pabe a , p a a = 2, p a b p' 
then (e° is shear and curl free) 

1 I, , „„ , 2 



c fa 



(63) 



D a eb = -Pab 
r 



--(Pabe c +Paceb), DVa 

r 



-e 
r 



(64) 
(65) 



so from (|48| ) 

D a D a (^0^) = -k 2 phys (R Ae O A t). 

To work out the l.h.s., we must first calculate D a Q£ ( |149[ ) and D a D a Qi (15C). These are then used 
to calculated : D a R Ae ( plD , D a B a R Af (jl5^ ), D a i?^ ([l53l) and D a D a O^~(p4D respectively. 
Putting these in (|150|) to find (p.56|) frow which we get : 



r 2 9r 



Or 



+ RA e 



1(1 + 1] 



(66) 



the spherical Bessel equation. Thus providing the PSTF derivation of the spherical bessel equation 
in terms of the irreducible representation. 

Now consider that we can choose any basis we like for the tensor basis here, independent of 
the spatial coordinates used. It is convenient to use the plane wave decomposition to get a parallel 
vector basis. We do this by writing 



R Ae (r) = J2R i A ) e =Y, R t(k,r)0 



(67) 



this expresses the tensor eigenfunction in terms of the monopole eigenfunction. When this is 
substituted into ( |156| ) we obtain the radial equation which has solutions that are spherical Bessel 
functions in the flat case : 



d_ 

Or 



,dRi 
dr 



+ 



k 2 r 2 

n phys 



Re = 0, to find R e (r) = m je(k phys i 



(68) 



where k 2 hys = k 2 /a 2 and ae are integration constants (the second set of constants for this second 
order equation vanish because we choose Ri(0) to be finite; the Neumann functions are not finite 
at r = 0). 



From (|6^) , fl67|) , and (p8|) we have found that the solutions to the Helmholtz equations give the 
eigenfunctions 



Qi(x" 



ae ^2h(kphysr)OfAO^ = a t ^2ji(k phys (k,r)r)L e (X). 



(69) 



k" 
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and we can set ai = (A^) _1 so that / QidQ = 1. 

It is important to notice that the functions L^{X) depend both on k a and on e a , and so for 
each k a is a function of (9,(f>), thus Q is indeed a function of all spatial coordinates. We can pick 
any direction k a to find the particular eigenfunctions 

Qf = (A^-^XWOO^O^. (70) 

associated with that direction. The general £-th eigenfunction is a sum of such eigenfunctions over 
a basis of directions k a qj: 

qk^) = E# ) - ( 71 ) 

k a 

Now we find Ae DA e Q ( |161| ) in terms of Ri (as show in the appendix using PSTF techniques, 



157-161) : 



M D WO = (-WK)°S^ ^ A *Q Al = (-k phva ) t Q i . (72) 
Putting this in the expansion ( [54] ) 

r^(x J ) = E^(*,fc)(-W)~' D <^>Q0O = EE^MX-Wr'D^Qn (73) 



gives the present version of fl52| ) 

directly analysing the coefficients in terms of these functions Q - given that we have 3-dimensions 
worth of variability so as to represent arbitrary spatial functions - with purely time-dependent co- 
efficients (parametrized by a vector k a ), 

4.1.3 Radial expansions and plane-waves 

Now consider the inversion]^] 



e (+ik a x a ) 



E^(r)Oj), <=► R Ae =Ai l Jdn {x) e+ ik "* a oW (75) 



by taking a Taylor expansion using k a x a = k p h ys re^ e h ^ we find 

6 _ a-1 y (+*)"( W r )" W /jnflWA 
n=0 

" ^ ^ n! ( n _^ + i)!( n + £ + i)H 47re ^^> ft 



n=0 



(76) 



on using (p7|) . Putting (|i20|) into (|7§) we find 



_ ( | .v A -i lT y (~i) m (W^ +2m (2m + 1)11 gW 



m _ n (2m + l)!(2m + 2^ + l) 
4ttI'+?VA~ 1 V (- 1 ) m (V^ r )^ +2m n (k) (77) 



14 Many treatments choose a particular direction for k: k a — 6% or similar and omit the summation. 
15 Contrast with (^2|): there the l.h.s is the spherical eigenfunction; here it is the plane one, expressed in terms of 
spherical ones. 
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This can be re-expressed as 

R Al = 4fK{+ifAj 1 j l {h phys r)0^} = (+iY(2£ + l)(3^j e ( k phys r)0^. (78) 

Hence 



Q flat (x\ e«) = e +lka *« = 4vr £ (+i) n j n (k p hy S r)OfyO^A~ 1 (79) 



n.fc" 



links the plane-waves to the spherical expansion This recovers, from Li(e a e' a ) = PePe(e a e' a ) the 
more usual 

e +ik phvs rX = £ i^o^ (x) = £(+^(2£ + l)j ( (Wr)P,(I) (80) 



i 

2 



4.1.4 Ir^l 2 for K = (flat) almost-FLRW models 

We now return to the relationship between the ta, ( and T£. Now from ( p-45|) , ( |148j ), (f79|) and ( |73[ ) : 

—-3 J dn k T £ (k,t)0% Y,i+ i ) n U^)Of~ ) 0% ) J~\2n + l). (81) 



n=0 



This can be reduced further ( 162 - 167]) , where upon using ( p.63|) one finds that 

, > r k' 2 rlk' 

r At = 4ttO { £ J Te (k',t)R t (k',x), (82) 

where \ = r /° an d by using R((k,x) = je(^ r ) we can identify this with ( |162j ), the equivalent of 
the multipole moments found using the explicit form of the plane-waves. 

In order to proceed further, from we construct the ensemble average ^] 



' = (4^o^>J ) /|^rKfc^^)i^K^^0/^rKA ; ,^)^(^r)\, 



47T 



(2vr)S 



(2tt) 3 ^ ' ' tv ' 'J (2tt) 3 
k 2 dk\r e (k,t)\ 2 . (83) 



What has happened here is that we imagine an ensemble of universes and we use an ensemble 
average rather than the space average; we have to do this because T£ is not square-integrable, i.e., 
we cannot use the r.m.s value as we cannot integrate the square over the space in general; this 
operation is not well defined. In order to deal with the ensemble average over the mode coefficients, 
i.e. { T e T m)i we assume the perturbations to be fairly homogeneously spread throughout the 
space and not confined in a particular region, and assume that there are no correlations between 
perturbations with different wavenumbers. Here the Gaussian assumption is useful, we have that 

(r £ (k,t)r e (k',t)) = (2irf5(k' - k)\r e (k,t)\ 2 . 



3 Remembering that 
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4.1.5 \t(\ 2 for almost-FLRW models 

It is useful to notice that an alternative although equivalent avenue of approach is also possible 

At , 



proceeding directly from the mode expansion in Gi[Q] = 0,AQA e ] this can be used for any K. 



Notice that as before at the observer at ccq, 

r(4,e a ) = Y^r Al O Al , r Ae (4) = Aj 1 f dOO Af r(4, e a ). (84) 

i J 

Now in the spatial section in general we can write 

r(x\e a {k) ) = E^O^), = AJ 1 [ d^O%r(x\e a {k) ). (85) 

e J 

Now the point Xq is chosen at some earlier time in a spatial section, with radial direction vector 
for FRW models we can consider e? \ and e a to be equivalent and write 

r^(4) = A7 1 / ^ (x) 0gr(4, ef x) ). (86) 
Remember that in some spatial section (where the integral relations ( |168 ,16£, 17C| ) are useful) 



r(x\e a ) = J2ri(k,t)0^ ) Q Ae . (87) 

Now from 



'(xr 

Lk v 



we identify O/M with O^, to find 



T A t =^2n(k,t)QA v (89) 



k a 



which means that for K = (it^ oc ji). 

jA e r Ae ) = Y,\^\ 2 (Q^Q Ae ) = ^mT,\^\ 2 - (9°) 



k" k 

What is particularly useful about the last method of calculation is that for constant K surfaces 
it can be shown on using relation ( |174| ) recursively that 

t 

(91) 



i 

/ dfl k (\Gt[Q]G m [Q]\) = 4xEffi , with S? = J] K)' 



and (| ... |) denotes the ensemble average of the eigenfunctions. What is useful here is to notice 
that (using the normalisation for K = 0), 
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This means that for K = 

J dQ {k) (|G$G*|) = ^^Af. (93) 

Now we can extend the above to FL models with K ^ 0. The mean-square for constant K is 
found by modifying the normalization after identifying the K = normalizations in (170) and ( B^ ) 

^Aj At ) = £ N 2 (Qa,Q A = (47r)/%£ |t,| 2 - 2 , (94) 

k v k 

or alternatively keeping the form of the mean-square in (^) by redefining the mode function 
expansion pj] by a wavelength-dependent coefficient; then 

T{x\e a ) =J2n(k,t)M e [Q], and M e [Q]=Ei(k,K)G e [Q], (95) 

defines the new coefficients fg(k,t) so that 

\f e \ 2 = \r e \ 2 Ej . (96) 



Because we have redefined the mode functions in (Pop, the form of the equations for K ^ is the 
same as in the case K = 0. However the coefficients are different because they are from a different 
expansion. Using the results from the K = case, either (|S2^) or ( |162| ), we have that 

(^r Ai ) = ^Pi J ^\ n (k,t)\\ (97) 

Hence, on using Ci = A _1 C£ and (rA e T Ae ) = Ci(2l + 1), we now have 

oo 

(r(x ,e a )r(x ,e' a )} = £(2£+ lrV^r^O^O'^ (98) 



which can now be written from (|97|) and ( pq ) as 

/ 1 \ 00 f°° dk 

(rK,e>K,e' a )) = J £(2^ + l)" 1 ^ ^ T ^\f e (k, t)\ 2 Q Ae O 



DO 



1=0 

thus reobtaining the results of White and Wilson : 

(3j f°° dk 



/ 1 \ r dk 



One needs to be careful here with the factors of (2£+l). Equation ( |99| ) relates the amount of power 
in a given wavenumber, \Tg(k, 770 ) | 2 , given the intersection of the null-cone fixed at the observer x , 
on a angular scale £ given that the angular correlations are found on the scale X = e a e' a , i.e. , X 
is the separation between measurements. 
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5 Conclusions 



We have given here a comprehensive survey of the CGI representation of CBR anisotropics in 
almost-FLRW universes, and related this formalism to the other major formalisms in use for this 
purpose at the present time. 

This paper has been concerned with algebraic relations: specifically the Multipole, (e.g.for ta ( ), 
and Mode, (e.g. for Tg), formalisms and the relationship between these. Where possible the multipole 
moments have been treated for a spacetime with generic inhomogeneity and anisotropy but small 
temperature anisotropies. The mode moments however are only meaningful in the restricted class 
of almost-FLRW universes. 

The subsequent papers in the series consider the differential relations satisfied by the quantities 
mentioned here [|14|] , and will show how both timelike and null integrations are used to lead to the 
standard results in the literature. Taken together, this will be is an ah initio demonstration of the 
way the different formalisms in use, and their major results, can be obtained from a single CGI 
approach, as well as providing the natural extension of the usual results into the non-linear (exact) 
theory. 
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A Spherical Harmonics 

A.l Basic relations 

A Spherical Harmonic (SH) Y^ m (9, (/)) is related to an Associated Legendre Polynomial (ALP) H, |TJ, 

Y im = c^e^P 1 " 1 (cos 9), (101) 

[(<-m)/2] 

= C im (e^sin0) m £ A lm \cos 9Y~ m ~ 2j V m > 0. (102) 
i=o 

Here, 

ctm _ W + W-my. and ^ _ J-VW-W {m 

v ' 4vr(£ + m)! ' 2 l j\{t - j)!(£ - m - 2j)\ v ' 

aloong with 

Y lm = (_i) m y'M* V m < 0. (104) 

Now we can relate the SH, Y m , to the direction vector product e Ae , 

Y* n = y%e A '(0,<l>), (105) 

where following H (from making the substitution e 1 + ie 2 = e 1 ^ sin and e 3 = cos into the above 
relation) 

[(t-m)/2] m £-2j 

j=0 k=0 p=m+l 

9=1 

Furthermore, it can then be shown that from 

m=+l 

f = Y J F Al e A ^ and F At = £ a^i™ (107) 

This is not unexpected. 

A. 2 Consequences 
A.2.1 Closure 

oo l=+m 



£ E Yi m (Sl')YUV>) = S(n - Of). (108) 

e=oe=-m 



A.2.2 Addition 

e=+m 



£ ^m(n')^m(n) = = a^o^o'^. (109) 

i=—m 
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A. 2. 3 Orthonormality 

/ dQY£i m i(Q')Yi m {Q) = 6 u >$mm>- (110) 

J4tt 

A. 2. 4 Matching plane-waves to spherical harmonics 

e (ikx) = A 7T J2 ie Mkx)Yi m Y , e ,- m = 47rJ2i e ji^7 1 Ae O' Ae , (111) 

i i 

i i 

B Multipole relations 

B. l Properties of e Ae 
B.l.l Normalization 

The normalization for e Ae is found from || , for odd and even £ respectively : 

_ f e A2 *+idU = 0, and —[ e A ^aTi = — !— h^ 2e) . (113) 
4W 4w 4vri 4 ^ 2£+l y J 

From which (contracting with h(A 2 e)) ^ can ^ e snown that 

^(^ = (21 + 1) (114) 
this can also be shown algebraically (7j. 

B.l. 2 Orthogonality 

From H we also have that 

/ e A *e B ™dn = — h( A * B ™\ (115) 

Jatt £ + m + 1 

if £ + m is even, and is zero otherwise (this follows from the above because e Ae e Bm = e Ae+m on 
relabeling indices: b\..b m — » a^ + i..O£ +m .) 

B.l. 3 Addition Theorem 

From (^) it follows that 

£=o f i cosp 

where X = cos (3. It may be useful to compare these to the relations for standard spherical 
harmonics, which are given in the appendix A. Note that 

/ dQe A "e' An = [ dQ( e V a )" = 2tt f ' dXX n = { ° 4?r V n °*J , (117) 

J4vr JA-K J -I [ ^+T V n eVen 

where the integral is taken over e a with e' a fixed. 



21 



B.1.4 Orthogonality of Al 

The orthogonality conditions can be found from 



£ J dn(F Ae O A *)(F Bm O B -), (118) 

see ||. Here F Ae are arbitary PSTF harmonic components of some f(e a ,x l ). Using flllEQ , (PI), 
and ([l^) we find 



J <mO A *0 Bm = 5 l m A e h^ {Be) with A, := ^ + ^ ^ . (119) 

where h^ At ' i B( \ = (61 " " " h ai \ e ) an d From this it follows that 

e B -h^\ Bl h Bn ^ = e^(+l) n - £ = A *. (120) 
It should also be noticed that from (|119|) , 

h^\ Ae) = (2£ + l), (121) 

can be also shown algebraically 0Qj 

Using these relations we obtain the inversion of the harmonic expansion (g): 

00 „ 
T(x i ,e a ) = Y,T Ae (x i )0 Ae O r A£ (x i ) = A7 1 / r(x\ e a ) . (122) 

B.1.5 Addition of Af 

The addition theorem for A( can be found from 

[1/2] [1/2] 

At O' Al = E E B ek B ek ,h^h {A2kl e A ^e' Ae) . (123) 

fc=0 fc'=0 



The resulting polynomial 



[1/2] 

L e (X) = A *0' Ae = £ B Em X e ~ 2m , (124) 

m=0 



is the natural polynomial that arises in the PSTF tensor approach. 

B.1.6 Double Integrals 

First, note that 



0' Ai e' Bn fdnO A 'e Bn = J r dnO' \0 At e> ' Bn e h ' 



4vr n\(n-£+l)\\ 
f e n + £ + l{n-£ + l)!(n + £ - 1)!! (125) 



Note the contrast with (114) 
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on integrating a Legendre polynomial, where n > / and (n — £) is even, so we can write n — £ = 2m 
for m an integer 0. Consequently, on remembering 

h^\ Bn h B ^ e) 0' Al e' B - = 0' Ae e' A * = fa, (e'V a )^ = +1 , (126) 

we find 

= ^"• ( „_^ ( 1 + + 1 |; 1) ,, ^»,b,^_. ) . (127) 

Here m are positive integers. Also we will need 

| dn k of t J dtivofy = Jdn k J cKiuO^ofy =4ttJ dn k o { £of k ? } 

= 4*6? Aeh^WSW -<%>)■ (128) 

C Mode relations 

C.l The curvature modified Helmholtz equation and the mode recursion rela- 
tion 

By successively applying the background 3-space Ricci identity, 

I K 

V abAe Q - D baAe Q = + E ^2 (^m " 5b aK nb ) V Ae Q, (129) 

n=l 

where = a\...a n -\b n a n+ x...ai, i.e. , the sequence of £ indices with the n-th one replaced with a 
contraction, the following useful relations are found: 

e^O A ^B At+1 Q = e ^O A 'B Ae Q-~^^0 A ^D Ae _ 1 Q, (130) 
h^O A ^D Ae+1 Q = (-k^ + l^^^jO^B^Q, (131) 



A ^ CJ CQ = (-k 2 phys + --£(£ + 3)) A *D Ae Q (132) 



A 'D cAe Q = A 'D AeC Q + --£(e-l)Oc e - 1 B A ^ 1 Q. (133) 
Now consider Ae B c cAe Q from ( ggg) we find 

G A ^ c cAt Q = O a *D\ c Q + —£(£ - lJO^DV.Q (134) 

where the first term on the left of the equality above, ( [L34| ), can be reduced to one in terms of 
Ae D A( Q using ( |132[) while the last term can also be rewritten in terms of At D A ,Q. Now on 
dropping the At and making the identification of Q A( = {—k p hy S )~ l B (A e )Q it is then shown that 
the Qal satisfy the curvature-modified Helmholtz equation 

B a B a Q {Ae) = (-k 2 phys + ^£(£ + 2))Q {At) , (135) 



18 This can also be seen from using F Al h {A * Bn) = (n!/(n - l)\\(n + i - iy.\)F Ae h Ae{Be h B "- e \ the definition of Ae 
in terms of e Af and evaluating J dD.F Ae O Ai e Bn for F AlL PSTF. 
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i.e., the Helmholtz equation with modified wavelength using as before k p h ys = k/a 



kj = (K£(£ + 2) - k 1 



(136) 



On using ( |135| ) and taking the PSTF part of the lower indices, D a D( a Q^, and using the PSTF 
tensor relation (|2l|) to find : 




(137) 



On substituting the first two relations ( |13C| , 131) into the recursion relation for the PSTF tensors 
|D, we find 

1)2 



e a B a [G e [Q}} = +k phys 



(2l + l)(2£-l) 



1 



K 



(r - 1) Ge-i[Q] - G i+1 [Q\ 



(138) 



C.2 Legendre tensors and the PSTF tensors 

One can immediately make the connection between this formulation and the one usually used in 
terms of Legendre tensors, and see that the Legendre tensors used in Wilson |l8}| in the coordinate 
basis (indicated by late romans) can be related to irreducible representation O e in terms of its 
associated tetrad frame Ef = |-u a ,e^j. The direction vectors e a in the triad with components 
are related to 7\ the direction cosines used in the Wilson-Silk coordinate basis treatment : 



P$ 2 -- Ae (f) = m- 1 el\eZ--<0 A t(e a ). 



(139) 



This connection to the Legendre polynomials can be seen by using the relation between spherical 
harmonics and the PSTF tensor along with the addition theorem for the PSTF tensors : 

Ae = yfaY lm (Q), (3 e Pi(e a e' a ) = A '0' Ae . (140) 

Multiply (H) by (3^ and use (p9|) and @ to find the recursion relation for the Wilson-Silk 
[17, [D| Legendre polynomials 



■pil—H 



{21 + 1) 



7 (#) 



^_ ^(ii»2 p*3— k+i) 



{£ + !) 1 ' W {I + 1) ' * V-V • 

It is seen that the recursion relations (gp for the irreducible representation O 1 can be reduced 
to that of the Legendre polynomial. This links the CGI-PSTF approach to the usual GI-Legendre 
tensor approach. 



D Plane-waves, spherical-waves and mode functions 
D.l Plane- wave and mode function relations 

We consider only flat, K = 0, universes at present. Each set of harmonic functions Q[k){ xa ) 
satisfying (48) which has associated with it kphys = k/a, the physical wavenumber, a variation 
vector field, q a , and a direction e a ( e a e a = 1, e a u a = ) determined by 

2 



D a Q = Qq a , q a = qe a , 



q q q a 



(142) 

19 The vector e a defined here is in general different from that associated with the angular harmonic expansion in 
(a). When ambiguity can arise, we explicitly put in the fc-dependence : g" fe j, to signify both this dependence and the 



definition of e a from (142) : thus strictly we should write, for example, D a Q(k) = Q(k) q ( a k) = Q(k)q {k) ei k) . We will 
suppress the k when this causes no ambiguity. 
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the first equality defining q a (x l ) (but not necessarily so as to factor out Q) and the second splitting 
it into its magnitude and direction. It follows that 

D a D a Q = Qq 2 + QD a q a = Q(q 2 + e a D a q + qD a e a ) (143) 

so that (|48| ) becomes 

q 2 + e a B a q + qD a e a = -k 2 phys , B a q a = -q 2 - k 2 phys . (144) 

(k) 

Using the K = plane -wave eigenfunctions with associated direction vector e a : 

Q( x \ e (fc))l/Zat = ex P {-ik p hyse { a ] x a } , (145) 

where k p h ys (k,t) = k/a, expresses the temperature anisotropy ([|) in terms of its plane-wave spatial 
Fourier Transform ([58]). In this case 



q a = -ik phys e[ k \ D a q b = = T> a k phys = D a ef \ q 2 = = -k 2 hys (146) 



holds in equation ( |142| ) and ( |143| ) respectively we find (from ( p.45[ )) 

^k))\flat = (-ik phys Y0f e Q{x\e {k; 



B {Ae) Q(x\e a {k) )\ flat = (-ik phys YO^Q(x\e a {k) )\ flat (147) 



(k) 

where the A[ are the PSTF tensors associated with the direction in the tangent spaces on 



the spatial section. Thus from (|49| ) 

{Q{x\e a {k) )\ flat ) {Al) = {-lfO^Q{x\e a {k) )\ flat (148) 

and (H) to find ©. 

D.2 Radial expansion and mode function relations 

D a Q e = D a (R Ai O A *) = Ae D a R Ae + R Ae D a O A *, (149) 

which implies that 

D a D a Q, = (D a D a R Ae )0 A * + 2(D a R Ae D a O A *) + R Ae (D a B a O A ^. (150) 

Now we need to work out (D a R Ae ), B a B a R Ae , B a O Ae and B a B a O Ae , say (a), (b), (c), (d) respec- 
tively. Calculating (a) : 

D^(r) = ^D a r = ^e ai (151) 
or or 



hence (b) follows: 



Next, (c) is: 



D a D a ^(r) = ^i + ^. (152) 



B a O Ae = - Pa ^O Ae ~^ e a D a O Ae = (153) 
r 
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which gives (d) : 

D a D a O Ae 



-(D a Pa ^)0 Ai -^ + -p a ^(D a O Ai -^) - ^D a (r)p a ^O Ae -^ 



-2£ 



l} + l{t - 1) 



(154) 



Now put these in ( |150| ) to find, 



D a D a Q e 



8 2 R Af , 2 0i2 A 



+ 



+ #A 



^ + 1) 



Af + 2 



A, 



kphysRA e O Ae , 



which simplifies to 



13/ 2 <9i?A f 



Af + R Ae 



l(£ + l) 



kphy S RA e O Al . 



(155) 



(156) 



Now we need to find 0,AD( A( }Q in terms of Rq, then relate the resulting rodrigues formulae of 

Re to get Ri in terms of Rq and hence Q A( in terms of Re. In this regard, now in (|67D, D a e^ = 
(see (p3D) which implies D a oJ ) = 0. Thus 



dRe{k,r) {k) _ dR (k,r) 



We can find Q(a £ ) from (49) obtaining 



<9r 



(157) 



D <^)Q = E D(^)^Og) = E D <A < ) J R m (^^)OgiO ( ^ ) mj = B {Ae) R (k,r), 



/ (PRo , 1 aRo 



D 



3 2 fi 

(Aj-3 I ^ r 3 e a*-2 e a*-i e a £ ) + ~ Q r 2 



2 5fip _l3flo m s \ 

g r P««-2a£-i e af> ^ r \ L> a e _ 2 r )Pa e _ 1 a e ) J ■ 



r dr 



(158) 



Now we note that 



1 d 

r 



(159) 
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which follows from p a6 e'< a e /6 > = h ab O' ab - e a e b O' ab = -O ab O' ab and that 

d_ 

r dr 



r e fid 

Ri = u— -z«z) (160) 



phys 

Used in (|159|) this gives that 



(x) 

E Mode mean square relations 
E.l Flat relations 

k 2 dk ' 



OfoD {Ai) Q = (-k phys ) e O A *oWR £ . (161) 



^ ' n=0 

" [ k 2 dkr £ (k,t)j £ (k X ). (162) 



2vr 2 

Equivalently from 

r(x\e a ) = Y,r m Qm, Qm = Y, R c m D m O^Of^ , and R Cm D m = Ri(k,r)h Cm D m . (163) 



(k) ~(x) 



Invert the multipole expansion 



t a , = A7 1 / dfi fc O« I £ r m Rc mBm Of k yO^ \ , (164) 

to find, on using (|12&j ) , that 

r Ae = ^Af^J^dkW^k^^Rc^O^^TAeh^}, (165) 



and hence that 



TAe = 2^J k ' 2dk ' T ^ k 'i O^QOg • (166) 

Using ( |163| ) this becomes 

TA t = ^0 { l] I k' 2 dk'r e (k', t)R e (k', X ), (167) 



E.2 Constant curvature relations 

We now have that 



J dn {x) 0^ } G e [Q] = 5 e m A e Q A *, (168) 

J dn (x) G e [Q]G m [Q) = 5TQ A "QaA^ (169) 
dn {k) Ge[Q]G' m [Q] = J dn^of^ReO^O^RmO^ = 8 t m A t lQO ( go'fa 

ldn ik) G e G m = 5i l A 2 ^±^-Rj. (170) 
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Furthermore, we have that from the recursion relations (| 

e^ x) D a [Gy = +k phys [ajG e ^ - G e+1 ] (171) 



where 



and using |17) 0, 



^ 2 = (2. + i^-i) ( 1 -f (£2 - 1) ) (172) 

| | dx"e a {x) D a [G e Ge-i] = (173) 
dn<*> (|G n | 2 ) = {a n f f dn {k) (|G n _!| 2 ) . (174) 



to find 



Here we have denned the mean square to pick out the power spectrum which is a function only of 
the absolute value of the wavelength for a Gaussian distribution (there is no directional dependence, 
the modulus is only dependent on the wave number). 

F 1+3 Ort ho- Normal Tetrad relations 

We use an orthonormal tetrad approach (cf. 0). Consider an orthonormal tetrad basis E a with 
components E^(x J ) relative to a coordinate basis; here indices a,b,c... , that is early letters, are 
used for the tetrad basis, while late letters i,j,k,... are used for the coordinate basis. The differential 
operators d a = E l a di are denned by the inverse basis components 

E?Ei = 4 O E?El = 5 a b . (175) 

The tetrad components of a vector X 1 are X a = EfX 1 , and similarly for any tensor. Tetrad indices 
are raised and lowered using the tetrad components of the metric 

g ab = gijElEi = diag(-l, +1, +1, +1) , g ab g bc = 5 a c (176) 

the form of these components being the necessary and sufficient condition that the tetrad basis 
vectors used are orthonormal, which we will always assume. 

For an observer with 4- velocity u a , there is a preferred family of orthonormal tetrads associated 
with u a i.e. a frame for which the time- like tetrad basis Eq is parallel to the velocity u a . In such 
a tetrad basis 

u a = 51 (177) 
hob = diag(0, +1, +1, +1) (178) 

All our work is based on such a tetrad, which leads to a preferred set of associated rotation 
coefficients. In paper 1, the form of these rotation coefficients is unimportant, so we defer their 
consideration to Paper 2. The issue for the present is that we have a preferred family of local 
orthonormal frames at each point (usually matter flow aligned), and carry out our algebraic analysis 
of observational quantities relative to that orthonormal frame. 



°The idea is to use this to fix the normalization of J dfl^)Gt[Q] = Di J df2(fc)Q(:r l , e^)0^0 



A, ■ 
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